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ABSTRACT 

We present optical photometry and spectroscopy of the X-ray transient XTE 
J1859+226, obtained during outburst and its subsequent decay to quiescence. Both 
the X-ray and optical properties are very similar to those of well-studied black hole 
soft X-ray transients. We have detected 3 minioutbursts, when XTE J1859+226 was 
approaching quiescence, as has been previously detected in the Soft X-Ray Tran- 
sients GRO J0422+32 and GRS 1009-45. By 24 Aug 2000 the system had reached 
quiescence with R=22.48±0.07. The estimated distance to the source is ~11 kpc. Pho- 
tometry taken during quiescence shows a sinusoidal modulation with a peak to peak 
amplitude of about 0.4 mag. A period analysis suggests that periods from 0.28 to 0.47 
days are equally possible at the 68% confidence level. The amplitude of the quiescent 
light curve and the relatively low ratio of X-ray to optical flux, indicates that the 
binary inclination should be high. The measured colours during the outburst allows 
us to obtain the basic properties of the disc, which agrees well with irradiated disc 
model predictions. 

Key words: black hole physics-binariesxlose-stars: individual: XTE J1859-I-226-X- 
raysistars 



1 INTRODUCTION 

Soft X-ray transients (SXTs) are a subclass of low-mass X- 
ray binaries (LMXBs) that are characterized by episodic X- 
ray outbursts (usually lasting for several months) , when the 
X-ra y luminosities can increase by as much as a factor of 
10^ ( |yan Paradijs fc McClintock 1995| ). The observed opti- 
cal flux is generated by X-ray reprocessing in the accretion 
disc and the companion star. These outbursts recur on a 
timescale of decades, but in the interim the SXTs are in a 

state rif qnipgpPTinp anH thp r.ptrr-a1 pmiggrnn ra HnminatpH 



best opportunity to analyze the properties of this star and 
obtain dynamical information which eventually enables us 
to constrain the nature of the compact object. There are 
currently 17 SXTs with identified optical counterparts, with 
13 dynamically studied black-holes and 4 confirmed neutron 
stars. 

The X-ray transient J1859-I-226, was discovered by the 
All-Sky Monitor {ASM) on board th e Rossi X-Ray Ti ming 
Explorer [RXTE) on 1999 October 9 (|Wood et al. 1999|) and 



subsequently identifie d with a R=15.1 variable star (dar- 



by the ladiation of the faint compauion star. This off-firs the navich fc Quinn 20001). Spectra during ou tburst showed typ 
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ical LMXB features ( [Wagner et al. 1999| ). The dereddened 
optical energy distribution was represented by a steep blue 
power law flattening in the ultraviolet which suggested a bi- 
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nary period <1 day. Although no persistent coherent mod- 
ulation was detected, a low amplitude (1 per cent modu- 



Table 1. Log of photometric observations (Oct 1999 -May 2000) 



latio n) with a period of 22-23 min was seen (Hynes et al 
1999). subsequent photome try revealed other m odulations 



with p eriods of 0.2806 days (Ucmura et al. 1999D and 0.3813 
days (Garnavich et al. 1999), but none of them could be 
confirmed. Nine months after the outburst, but still not in 



quiescence, a new 0.78 day period was reported (McClintock 
et al 



200C) 



rowed a 0.2 



Light curves folded on this period sh 
mag primary minimum (interpreted as a partial eclipse of 
the accretion disc by the secondary star) , and a shallow sec- 
ondary minimum (explained as an eclipse of the star by the 
disc) . 

Here we present our extensive optical photometry of 
the entire outburst, the subsequent minioutburst and the 
first quiescent photometry. Some of these results have re- 
ceived preliminary announcements in Sanchez-Fernandez et 
al. 2000, Zurita et al. 2000 and Charles et al. 2000. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Photometry 

Our long-term monitoring campaign was carried out during 
the period Oct. 1999- Sep. 2000 with the 80 cm IAC80 and 
Im Optical Ground Station (OGS) at the Observatorio del 
Teide; Im JKT, 2.5m NOT, 3.5m TNG and 4.5m WHT at 
the Observatorio del Roque de los Muchachos; 1.5m OAN 
and 2.2m CAHA at Observatorio de Calar Alto and 1.52m 
telescope at Loiano. We obtained CCD images mainly in R— 
band but also some B,V, I colours. Integration times ranged 
from 30 s to 40 min , depending on the telescope, atmo- 
spheric conditions and the star brightness. The observing 
log with full details is given in Table |^. All the images were 
de-biased and flat-fielded in the standard way using iraf. 

We applied aperture photometry to our object and sev- 
eral nearby comparison stars within the field of view, us- 
ing IRAF. We selected several comparison stars which were 
checked for variability during each night, and over the en- 
tire data set. Calibration of the dat a was perform ed using 
17 standard stars from several fields (Landolt 200C), to con- 



struct a colour dependent calibration. Frames taken during 
quiescence with the WHT and good seeing conditions, re- 
vealed the presence of a faint nearby star ~1.4 arsec North 
of the target (see Figure [l|). We applied psf photometry using 
the IRAF routine DAOPHOT ( [Stetson 198"^ ). We found a 
magnitude of 7?=23.05±0.04 for the contaminating star. We 
also recalibrated a set of 5 faint comparison stars from the 
previous calibrated set, using these images (see Table ^). 



2.2 Spectroscopy 

Spectroscopic observations of XTE J 1859-1-226 were carried 
out at primary outburst and during one of the miniout- 
bursts. The observing log for these observations is presented 
in Table |^ The first series of optical spectra of the source 
were obtained on Oct. 28-29, 1999, when 31 exposures were 
taken in the Isaac Newton Telescope (INT) at the Observa- 
torio del Roque de los Muchachos, using the EEVlO camera 



Date 




HiJjP / r LLLtil 


Te Igs cope 


OctlS 


7 


IxR 


0.8m IAC80 


Octl7 


9 


10xI,21xR,2xV,2xB 


0.8m IAC80 


OctlS 


10 


lllxR,lxV,lxB 


0.8m IAC80 


Octl9 


11 


106xR,lxV,lxB 


0.8m IAC80 


Oct20 


12 


1x1,1 10xR,lxV ,lxB 


0.8m IAC80 


Oct21 


13 


2xR,lxV,lxB 


0.8m IAC80 


Oct23 


15 


IxR 


0.8m IAC80 


Oct28 


20 


3xR 


0.8m IAC80 


Oct29 


21 


lxI,lllxR,lxV,lxB 


0.8m IAC80 


NovOl 


29 


128xR 


1.5m OAN 


Nov27 


50 


30xR,lxV,lxB 


0.8m IAC80 


Nov28 


51 


40xR,lxV,lxB 


0.8m IAC80 


DiclO 


63 


2xR 


0.8m IAC80 


Dicl2 


65 


4xR,lxV,lxB 


0.8m IAC80 


Dicl6 


69 


IxR 


0.8m IAC80 


Jan03 


87 


3xR 


2.2m CAHA 


JanOJf 


88 


3xV 


2.2m CAHA 


Jan05 


89 


IxB 


2.2m CAHA 


FebOl 


116 


IxR 


0.8m IAC80 


Feb04 


119 


lxI,lxR,lxV,lxB 


Im OGS 


Feb05 


120 


lxI,lxR,lxV,lxB 


Im OGS 


Feb06 


121 


lxI,2xR,lxV,lxB 


Im OGS 


FebOl 


122 


lxI,lxR,lxV,lxB 


Im OGS 


FebOS 


123 


lxI,lxR,lxV,lxB 


Im OGS 


Feb09 


124 


lxI,lxR,lxV,lxB 


Im OGS 


FeblO 


125 


lxI,lxR,lxV,lxB 


Im OGS 


Febll 


126 


2xI,2xR,2xV,2xB 


Im OGS 


Febl2 


127 


lxI,2xR,lxV,lxB 


Im OGS 


Febl3 


128 


2xI,2xR,lxV,lxB 


Im OGS 


Febl4 


129 


3xI,2xR,lxV,lxB 


Im OGS 


Feb 15 


130 


2xI,lxR,2xV,2xB 


Im OGS 


Feb20 


135 


IxR 


Im OGS 


Feb27 


142 


IxR 


0.8m IAC80 


Mar03 


145 


IxR 


0.8m IAC80 


Marie 


160 


IxV.lxR 


0.8m IAC80 


Mar23 


167 


IxV.lxR 


0.8m IAC80 


Mar24 


167-8 


lOxR 


Lowell 






lxV,lxR 


0.8m IAC80 


Mar25 


168 


IxR 


Lowell 


Mar27 


170-1 


30xR 


Lowell 


Apr03 


178 


lxV,lxR 


0.8m IAC80 


Apr04 


179 


lxV,lxR 


0.8m IAC80 


Apr05 


180 


IxR 


0.8m IAC80 


Apr 10 


185 


IxV.lxR 


0.8m IAC80 


Apr 16 


191 


IxR 


0.8m IAC80 


Apr25 


200 


IxV.lxR 


0.8m IAC80 


Apr28 


203 


2xR 


0.8m IAC80 


May 16 


221 


2xR 


0.8m IAC80 


May 17 


222 


IxR 


0.8m IAC80 


May23 


228 


2xR 


1.52 Loiano 


May24 


229 


2xR 


1.52 Loiano 


May25 


230 


2xR 


0.8m IAC80 






3xR 


1.52 Loiano 
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Table 

2000) 



1 - continued Log of photometric observations (Jun-Nov 



Table 2. Log of spectroscopic observations 



Date 


HJD^*'> 


Exp/ Filter 


Telescope 


JunOl 


237 


37xR 


Im OGS 


Jun03 


239 


34xR 


Im OGS 


Jun04 


240 


64xR 


Im OGS 


Jun06 


242 


69xR 


Im OGS 


JunOl 


243 


28xR 


Im OGS 


Junll 


246 


IxR 


0.8m IAC80 


Junl2 


248 


IxI.lxV 


2.5 NOT 


Junl4 


249 


IxV.lxR 


0.8m IAC80 


Jun24 


260 


IxR 


0.8m IAC80 






18xR 


Im OGS 


Jun25 


261 


17xR 


Im OGS 


Jun28 


264 


17xR 


Im OGS 


JulOl 


267 


lxV,lxR 


0.8m IAC80 


Jul02 


268 


2xR 


1.52 Loiano 


JulOS 


269 


IxR 


1.52 Loiano 


JuW4 


270 


lxV,lxR 


0.8m IAC80 


JulOS 


271 


336xB 


2.5 NOT 


JuW6 


272 


lxV,lxR 


0.8m IAC80 






887xB 


2.5 NOT 


Julll 


277 


IxV.lxR 


0.8m IAC80 


Jull4 


280 


lxV,lxR 


0.8m IAC80 


JullS 


281 


lxV,18xR 


0.8m IAC80 


Julie 


282 


IxR 


Im JKT 


Jull7 


283 


IxR 


0.8m IAC80 


Jul 19 


285 


IxR 


0.8m IAC80 


Jul20 


286 


IxR 


Im JKT 


Jul21 


287 


IxR 


0.8m IAC80 






IxR 


1.52 Loiano 


Jul22 


288 


IxR 


0.8m IAC80 


Jul23 


289 


IxR 


1.52 Loiano 


Jul24 


290 


IxR 


1.52 Loiano 


Jul27 


293 


IxR 


1.52 Loiano 


Jul30 


296 


168xR 


Im OGS 


Jul31 


297 


59xR 


Im OGS 


AugOS 


302 


IxR 


0.8m IAC80 


AugOe 


303 


IxR 


0.8m IAC80 






IxR 


1.52 Loiano 


AugOl 


304 


IxR 


1.52 Loiano 


AugOS 


305 


IxV.lxR 


0.8m IAC80 


Aug09 


306 


IxV.lxR 


0.8m IAC80 


AuglO 


307 


IxR 


1.52 Loiano 


Augll 


308 


IxR 


1.52 Loiano 


AuglS 


315 


IxR 


1.52 Loiano 


Augig 


316 


IxR 


1.52 Loiano 


Aug20 


317 


IxR 


1.52 Loiano 


Aug21 


318 


IxR 


1.52 Loiano 


Aug24 


321 


IxR 


0.8m IAC80 


Aug25 


322 


IxR 


0.8m IAC80 


Aug27 


324 


IxR 


0.8m IAC80 


Aug28 


325 


IxR 


0.8m IAC80 


Aug29 


326 


IxR 


0.8m IAC80 


Aug04 


332 


IxR 


0.8m IAC80 


AugOS 


333 


IxR 


0.8m IAC80 


Sep27 


355 


6xR 


4.2m WHT 


Sep28 


356 


16xR 


4.2m WHT 


Octll 


374 


IxV 


2.5m NOT 


Nov04 


393 


6xR 


2.5m INT 


NovOS 


394 


5xR 


2.5m INT 



Date 


HJD('^ 


Number of spectra 


Telescope 


99 Oct28 


20 


1 


1.52m G.D. Cassini 


99 Oct28 


20 


7 


2.5m INT 


99 Oct29 


21 


24 


2.5m INT 


00 Julll 


277 


12 


4.5m WHT 



*HJD-2451460 




'HJD-2451460 



Figure 1. R band 1800 s image of XTE J1859+226. The field of 
view is 1.4 x 1.0 arcmin. XTE J1859+226 is indicated by an arrow 
and the contaminanting star by the letter C. The magnitudes of 
stars 1—5 are hsted in Table ^ 



on the Intermediate Dispersion Spectrograph (ids) (range 
3500-5000 A, spectral resolution 1.02 A/pixel). The expo- 
sure times ranged from 300-600 s for individual spectra. A 
single 2400 s spectrum was obtained with the 152cm G.D. 
Cassini Telescope at Loiano Observatory using bfosc. The 
slit width was 2.5" and the spectral range 3500-9000 A. The 
second series of spectroscopic observations were carried out 
on July 11 2001, during the second optical miniout burst de- 
scribed in section ^ We observed XTE J1859-I-226 using 
the ISIS spectrograph on the 4.2m WHT in Observatorio del 
Roque de los Muchachos, using typical exposure times of 
1800 s. Standard iraf procedures were used to de-bias the 
images and to remove the small scale CCD sensitivity vari- 
ations. One dimensional spectra were then extracted from 
the processed images using the optimal extraction method 
(Horne 1986). Wavelength calibration was interpolated be- 
tween contemporaneous exposures of a copper-argon arc 
lamp. No flux calibration was performed on the spectra, in- 
stead continuum normalization was applied. 
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OctOB Novl7 Decl7 Feb05 Marl6 Apr26 Jun05 Jull4 Aug24 Sep03 




40 80 120 160 200 240 280 320 360 

HJD-2451460 




640 



HJD-2448895 

Figure 2. Top: Temporal evolution of J1859+226 plotted as 'X-ray magnitudes' [24-2.5/og(XTE/ASM Count Rate)] and R-band 
magnitudes averaged for each day. Note that 1 Crab equals an xte/asm count rate of 75 counts/sec. A, B and C mark the epochs studied 
for photometric and spectrometric variability . The arrows indicate the dates when spectra were taken. Bottom: R band light curve of 
XTE J0422-I-32 plotted for comparison. The J0422+32 data were provided by E. Kuulkers. 
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Days 

Figure 3. Top: The decline of tlie main outburst of J1859+226 showing the subsequent minioutbursts. The three minioutbursts are 
marked as Jun 00, Jul 00 and Aug 00. Bottom: The minoutburst of J0422+32 plotted for comparison (asterisks). We have superimposed 
the J1859+226 minioutbursts on the same scale (filled circles). For both systems, the last point marks the day when quiescence is reached. 
The J0422+32 data were provided by E. Kuulkers. 



Table 3. R magnitudes for the contaminating star and 5 faint 
comparison in the field of J1859+226 



c 


1 


2 


3 




5 


23.05 
±0.04 


20.31 
±0.06 


19.41 
±0.05 


18.86 
±0.05 


20.15 
±0.06 


19.58 
±0.05 



3 LONG TERM BEHAVIOUR 

Both the X-ray and optical properties of XTE J1859+226 
were very similar to those of previously well-studied black 
hole SXTs, and so we embarked on a campaign of systematic 
monitoring of the decay Ught curve. Optical outburst Ught 
curves of SXTs tend to be fragmentary and generally cover 
only the main outburst. We note that the only SXT that 



has been extensively covered throughout the entire outburst 
and subsequent decay activity is GRO J0422+32 (Callanan 
et al. 1995, Chevalier & Ilovaisky 1995) . Figure ^ therefore 
compares the overall 1-yr Ught curve of GRO J0422+32 with 
of XTE J1859+226 in optical (i?-band) and X-rays (2-12 
keV), since the X-ray turn-on of 9 October 1999. The ba- 
sic properties of the outburst light-curve, are summarized 
in Table ^. They are similar to the Ught curves reported 
in Chen et al. 1997. The optical hght curve can be classi- 
fied as FRED: 'fast-rise exponential-decay', characterized 
by a smooth decline of 0.017 mag/day. At about the time 
that the X-rays dropped below the RXTE sensitivity limit, 
the optical light curve began a precipitous fall from i?~18.5 
to ii~22.5 (day 210 after outburst). About 12 days after 
the peak of the outburst, the X-ray intensity reaches a sec- 
ondary maximum which is not visible in our optical data. 
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Table 4. Basic properties of the outburst. 



OctOO Novl? Dec17 FebO.-) MarlB Apr26 JunOS Jul H AueS'1 Sep03 



Date of the outburst: 



Oct.1999 



X-ray light curve 



Morphology type: 

Instrument and E band: 

Flux at the peak: 

Rise timescale (rr): 

Decay timescale (t^): 

Duration of the rising phase (Tr): 

Duration of the decay phase (T^): 



FRED 

ASM(RXTE)/2-12 KeV 
0.25 Crab 
2.5 days 
34 days 
13 days 
155 days 



Optical light curve 



Morphology type: 

Peak magnitude 

Quiescent magitude 

Outburst amplitude (mqui^t-i^peak) 

Decay timescale (t^): 

Duration of the decay phase (T^): 



Exponencial decay 
No rise phase data 
R=15.1,V=15.3 mags 
R=22.48, V=23.29 mags 
AR=8.4, AV=7.8 mags 
103 days 
220 days 



A 'glitch' of Ai?~0.5 mag was detected in the expo- 
nential decay phase, at about Mar 25 (~165 days after 
the peak), and in the X-ray light curve about 5 days ear- 
lier (~160 days from the outburst). Similar glitches have 
also been observed in the optical light curve of A0620-00 
(Tsunemi et al. 1977), with a magnitude increase of Ai3~0.5 



mag about 150 days after the onse t of the outburst, and also 
in J04224-32 with AT/~0.3 mag (ICallanan et al. 19951). In 



X-rays, glitches are seen in the li ght curves of J04 22-I-32 



ight 

A0620-00, 1543-47, and 1124-683 jChen et al. 19971 ) 



Shortly after this rapid decline began, we have detected 
3 minioutbursts or small amplitude (relative to the primary 
peak) events superposed on a 'normal' decay profile, when 
XTE J1859-226 was approaching quiescence (see Figure. |^. 
The first one (Jun 00) occurs ~240 days after the peak, has 
an amplitude of 1.3 mag in R (with respect to the base level 
of Jun 2000) and lasts for ~20 days. After the minioutburst, 
the decay follows roughly the same decline as before. The 
second minioutburst (Jul 00) occurs ^265 days after the 
peak, reaches a maximum of i?~18 (AR~4.7mag) and lasts 
for ^^25 days. And the third (Aug 00) reaches a peak only 
~2 mags above the previous minimum. The source then de- 
cays linearly to R ~23 about 30 days after the onset on this 
last event. Minioutb ursts have only been previous ly detected 
in GRO J0422-t-32 (|Chevalier fc Ilovaisky 1995| ) and GRS 
1009-45 ( Bailyn 1992| ) but this may be a consequence of 
inadequate monitoring (or lack of sensitivity) once the out- 
burst is over. In JQ422-f-32, two principal minioutbursts were 
observed (see Figure j^. Both had an amplitude of AF~5 
mag and lasted for ~20-40 days. Here the minioutbursts 
reached the same level as the extrapolated light curve be- 
fore the precipitous fall and are followed by events of smaller 
amplitude, as occur in J 1859-1-226 Jul 00 minioutburst. De- 
spite these similarities, the shape of the long term light curve 
in J0422-I-32 is completely different. Here the rate of decay 
in the i?-band i s only 0.0056 mag/da y, which is a remark- 



ably slow rate (Callanan et al. 199E). We note that both 



the duration of the main outburst and the minioutbursts is 
approximately the same. 
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Figure 4. Long term variations of the V /, V -R and 

B V colour indices of J1859+226 are plotted here against the 

overall light curve (bottom). 



Figure g presents colour information of J1859-I-226 as 
a function of time. The mean V-R colour during outburst 
was measured to be ~0.27 and remains constant during 
the main outburst. The source does not show detectable 
colour changes during any given night. It is clear that the 
system reddens as the system approaches quiescence and 
the secondary's contribution increases. During quiescence 
V--R=0.81±0.11, which is consistent with the colour of a 
G9-K5 main sequence star. 

From 24 Aug the system has reached quiescence with 
a mean magnitude of i?=22.48±0.07 and y=23.29±0.09 
(measured on 17 Oct). This yields a total amplitude for 
the optical outburst of AV=7.8 mag, com parable to that 
observed for other SXTs (|Chen et al. 1997|). 



4 ANALYSIS OF PHOTOMETRIC 
VARIABILITY 

Photometric variability has been explored in our best sam- 
pled data set. This includes a group of 6 nights at the peak of 
outburst (18-20, 29 October 1999, 8 and 28 November 1999), 
2 nights during the maximum of the July minioutburst (5- 
6 July 2000) and 4 nights in quiescence (27-28 September 
2000 and 4-5 November 2000). Hereafter, we will refer to 
data taken during the main outburst, the Jul 00 miniout- 
burst and quiescence as Epochs A, B and C respectively 
and they are marked in Figure. |^. 
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4.2 Epoch B: 22 min QPOs 



K 15.6 - 



18/10/1999 


29/10/1999 
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1 1 1 1 

20/10/1999 


28/11/1999 



Figure 5. Optical light curves for XTE J1859+226 on different 
nights in Epoch A. Note the remarkable changes in the light curve 
shape during each night. 



4.1 Epoch A: Irregular Variability 

Representative light curves of XTE J1859+226, at the be- 
ginning of the outburst, are presented in Figure We note 
the remarkable changes in shape and amplitude with no 
evidence for repeatability. Intermittent modulations of this 
type during decline from outburst on timescales at or near 
the orbit al period have been s een in GS 200 0+25 (Charles et 
al. 199l[ ), Nova Muscae 1991 ( [Bailyn 199^ and Nova Persei 
1992 (Callanan et al. 1995; Chevaher and Ilovaisky 1994). 



Fourier analysis of these October data, shows three 
main peaks of comparable significance at 0.1600, 0.1907 and 
0.2353 days in the power spectrum. However, no clear modu- 
lation is obtained after folding the data on these periods. We 
also note that none of these periods have been obs erved by 
other auth ors. For example periods of 0.2806 days ( Uemura 
et al. 1999) and 0.3813 days ( [Garnavich et al. 199£| ) were re- 
ported during the first phase of the outburst. Nonetheless, 
the determination of orbital periods of X-ray novae during 
outburst has proved notoriously unreliable, e.g. modulation 
in the outburst light curve at a period much shorter than 
the orbital period has been seen in V404 Cyg (i.e. Wagner et 
al. 1990, Gotthelf et al. 1991). On the other hand it has been 
suggested that J1859-I-226, has a binary period shorter than 
1 day, since the dereddened optical spectral energy distribu- 
tion of J1859+226 during outburst can be represented by a 
steep blue power law fiattening in the ultraviolet, and this 
r esembles energy d istributions of other short-period SXTs 
( iHynes et al. 1999| ) . 



High time resolution (30 to 60 s) B-band light curves were 
obtained on 5 and 6 July when J 1859-1-226 was at the max- 
imum of a minioutburst (see Figure. |^. The light curve on 
Jul 5 shows an extended 'dip' of about 0.3 mag lasting for 
3 hr. During this dip, the source exhibited remarkable QPO 
flaring activity with ~20 per cent amplitude. On July 6, the 
QPOs has disappeared and instead the source exhibited a 
smooth sinusoidal modulation with ~0.10 mag amplitude. 
Studying the rapid variability seen on July 5, we find a pe- 
riod of 21.7±0.6 min using a periodogram analysis. This 
QPO timescale is very close to that reported by Hynes et 
al. 1999. The flaring ceased when the source returned to its 
predip level. 

A periodicity of 0.7 8 days has also been r eported 
from contemporary data ( McClintock et al. 200C), which 
is claimed to be the orbital period. Their light curve shows 
evidence for a deep primary minimum (interpreted as the 
eclipse of the disc by the companion star) and a shallow 
secondary minimum. Assuming that the 0.3 mag dip in the 
July 05 light curve is the primary minimum, it occurs at 
HJD 2451731.51, which is exactly 1.3 days before McClin- 
tocks's et al. primary minimum. Therefore the true orbital 
period would be a submultiple of 1.3 days, which au tomati- 
cally rules out the 0.78 d period (Zurita et al. 2000). 



4.3 Epoch C: Quiescent ellipsoidal modulation 

Data taken in quiescence during 2000 September 27 and 28, 
on the 4.2 m WHT, and 2000 November 4 and 5 on the 2.5 m 
INT, exhibit a sinusoidal modulation, with a peak-to-peak 
amplitude of about 0.4 mag. We applied a Chi-squared min- 
imization method in order to characterize the periodicity 
present in the data. Periodogram is shown in Figure |^ The 
periodogram shows several peaks although formally the peak 
with the highest signiflcance is found at 0.1594±0.003 days. 
Assuming that this is actually a double-humped ellipsoidal 
light curve, as is typical of quiescent soft X-ray transients, 
the orbital period would then be 0.3188±0.003 days (~7.7 
hours) . We note that this period is roughly a submultiple of 
1.3 days. This suggests that the variability (and in particular 
the 0.35 mag drop observed in the July minioutburst), might 
be related to the orbital period. However, periods from 0.14 
to 0.23 days (Port ~6.6 to 11.2 hrs ) are equally possible 
at the 68 percent confidence level. Clearly more and higher 
quality data are needed to determine conclusively which pe- 
riod is the true orbital period. 

We note that Filippenko & Chornock (2001) report a 
9.16±0.08 hr (0.382±0.03 days) modulation in the radial 
velocities of only 10 spectra, spread over two nights, ob- 
tained during quiescence. This period is consistent with 
our data and previous reports by Garnavich et al. (2001) 
and Sanchez-Fernandez et al. (2000). Since all the dataset 
are affected by aliasing, it is clear that more spectro- 
scopic/photometric observations are required to determine 
the true orbital period. 

In Figure ^ we show the quiescent data folded on the 
0.319 days period (top panel) and on the 0.382±0.03 days 
Filippenko's spectroscopic period (bottom panel). By com- 
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Figure 6. B-band light during Epoch B (JulOO minioutburs) of J1859+226 on 5 (top) and ' 
obtained with the 2.5m NOT. Note the flares in the Jul 5 light curve. 



Jul (bottom). The data points have been 



paring the amplitude of the curve with other SXTs in quies- 
cence, we conclude that the the binary inclination must be 
high. The different depths in the minima also support this 
idea since this is only evident when the inclination angle is 
high. 

An alternative way to estimate the orbital period is 
through the V amplitude of the outburst light curve (AV^), 
by the use of the empirical relation 

AV = 14.36 - 7.6 log Port (/ir) 



( Phahbaz fc Kuulkers 19981 ). Taking V(peak)=15.5 and 
V(quiescence)=23.3, we obtain Porb ~7.30 hr (0.30 days) 
which is consistent with the previous values. 



5 ANALYSIS OF SPECTRAL VARIABILITY 

Spectral variability has been explored in epoch A (Oct. 20- 
28) and in the epoch B minioutburst (Jul. 11). Days when 
spectra were taken are marked by arrows in Figure. ^ and 
are plotted in Figure. H. 



5.1 Epoch A 

The bluer spectrum displayed typical emission lines of X- 
ray transients in outburst: He ll 4686 A, the Bowen blend at 
4630-4640 Aand 3097- 4103 A, and Balmer emission from 
H/3 to H(5. The H/3 emission line is embedded in a broad ab- 
sorption, as is often observed in the Balmer emission lines 
from SXTs (e.g. Casares et al. 1995). The main interstel- 
lar features were the Can interstellar bands at 3933.66 and 
3968.47 A, (probably with some contribution from He ab- 
sorption at 3970 A), and the difluse interstellar band at 4430 

A. 

The only prominent feature in the redder spectrum is 
Ha emission. The main interstellar features were the blend 
at 5778, 5780 and 5797 A, the blend due to the interstel- 
lar Na I D lines at 5889 and 5895 A, and the interstellar 
band at 6280 A. The Hell and Balmer lines appear to be 
double peaked with separations of 300-500 km s~^, notably 
smaller than those detected in the outburst spectra of GRO 
J0432-I-22 (Callanan et al. 1995) or in the quiescent spectra 
of Nova Muscae 1991 or A0620-00 (|Orosz et al. 1994|). Dou- 
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Porb ~0 
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7. minimization spectrum of the quiescent data. The 
minima are marked, corresponding to Port ~0.32 and 
.38 days. 
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Figure 8. Optical photometry of XTE J1859+226 during the qui- 
escent state (Epoch C) folded on the 0.319 d period (top panel) 
and on the 0.382 d period found by Filippcnko & Chornock (bot- 
tom panel) and a superimposed sinusoidal fits. We have repre- 
sented points taken on different days/sites with different sym- 
bols: asterisks (27/09/2000 - WHT), circles (28/09/2000- WHT), 
squares (04/11/2000 - INT) and triangles (05/11/2000 - INT). 



ble peaks were not detected in the Bowen blend lines in our 
spectra, most probably due to the blending phenomena than 
to the absence of such structure. However, the Bowen blend 
could also arise from other emitting region in the system. 



5.2 Epoch B 

The second series of spectra of the source were obtained on 
July 11 2000, near the maximum of that minioutburst. The 
bluer spectra displayed prominent Balmer emission from H/3 
to H9. Bowen emission was not detected in these observa- 
tions. On the other hand, Cll emission at 4267 Awas present. 
The main interstellar features were again the interstellar ab- 
sorption at 3968 A, which affects the He emission line, the 
diffuse Call bands and the diffuse interstellar band at 4430 
A. The Balmer lines also appear to be double peaked in min- 
ioutburst, with typical separations of 500-700 km s^^, no- 
tably larger than the separations observed during main out- 
burst. The red spectra display Ha emission and Hel emission 
lines at 5877 and 6680 A. The He I line at 5877 is affected by 
the interstellar blend of the Na 1 D. The interstellar bands 
at ~ 5780 Aand 6280 Aare also detected. 

The only observations of previous miniout burst spec- 
tra of SXTs are those from GRO J0422-I-32 (jCasares et| 



al. 1995; Callanan et al. 1995) and GRS 1009-45lBaityn 



fc Orosz 1995). The spectra from GRO J0422-f32 showed 



broad shallow H/3 absorption, and Ha evolving from ab- 
sorption to emission on a timescale of 3 days. The spectra 
from GRS 1009-45 showed weak Hq and weak H/3 absorp- 
tion. In none of them evidence of He I is found. Casares et 
al. 1995 also obtained phase resolved Ha and H/3 spectra of 
J0422+32 during its 1993 minioutburst. Balmer lines were 
embedded in broad absorptions whereas Hell was purely in 
emission and showed evidence of a large amplitude (~755 
km s^^) S-wave component. 



5.2.1 Variability of the emission lines 

To examine the variability of the emission lines during the 
minioutburst, the equivalent widths were measured in the 
individual spectra as a function of 0.319 d period (see Fig- 
ure. |l^. The equivalent widths of the lines were measured 
after normalization of each individual spectrum by its con- 
tiimum so that the variations in the lines are separated from 
variations in the continuum. Unfortunately, this study was 
not possible with the main outburst data, due to the poor 
S/N ratio of the individual spectra. Sinusoidal fits have been 
included in the plots as a reference. We can see that the Ha 
equivalent width is dependent on orbital phase. The same 
relation is less clear in H/3, probably due to the smaller am- 
plitude of the variations, but it is still present. The variation 
in the equivalent widths for both of them can be roughly 
reproduced with a sinusoidal function equal to the orbital 
period of the system, thus pointing to the changing visibility 
of a region of enhanced emission or bright spot, perhaps the 
splash point where the gas stream and accretion disc meet. 
Unfortunately, the uncertainty in our To prevents to define 
the absolute phasing in the EW curves. No clear modula- 
tion is seen in the EWs of He and Hell which have means of 
5.95±0.95 Aand 5.41±1.40 Arespectively. 
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Figure 9. Spectroscopic observations of XTE J1859+226 during Epoch A (main outburst)— figures a and b- and during Epoch B (July 
00 minioutburst)-figures c and d. 



6 DISCUSSION 

6.1 Distance estimate 



In the context of King & Ritter's model ( 1998), the expo- 
nential decay of the X-ray hght curve indicates that irra- 
diation is strong enough to ionize the entire accretion disc. 
Also, a secondary maximum is expected one irradiated-state 
viscous time after the onset of the outburst and this can be 
used to calibrate the peak X-ray luminosity and hence the 
distance to the source dkpc through 

d,p. = 4.3 X W-'t'/W^'f'^'F-'^\;'^' 

where Fp is the peak X-ray flux, the time of the secondary 
maximum after the peak of the outburst in days, Td the e- 
folding time of the decay in days, rj the radiation efficiency 
parameter and / the ratio of the disc mass at the start of 
the o utburst to the maximum possible mass ( ^hahbaz et al 



mated from the XTE count rate (250 mCrab in the energy 
range 2-10 keV) which corresponds to 1.7 x 10~* erg cm~^ 
s~^. Assuming 77=0.15 and /=0.8 we find dkpc = H- 



1998). In our case, rd=34 d, ts 



Alternatively, we can estimate the distance to the source 
by comparing the quiescent magnitude with the absolute 
magnitude of a main sequence star which fits within the 
Roche lobe of a 7.65 hr orbit. Combining Paczynski's (1971) 
expression for the averaged radius of a Roche lobe with Ke- 
pler's Third Law we obtain the well-known relationship be- 
tween the secondary's mean density and the orbital period: 
p = 110/P^j (g cm"''). Substituting for the orbital period 
of J1859+226 we obtain p=1.87 g cm~^ which corresponds 
to a KOV-KIV secondary star with absolute magnitude 
AIr ~6.1. The dereddened quiescent magnitude is V=23.27 
(using Av=1.80±0.07 as derived from the NalD line; see 
Hynes et al. 1999) which yields dkpc = H- This value per- 
fectly agrees with that obtained previously, although strictly 
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Table 5. Relative contribution of the disc in the R and V band 
(bold), assuming different spectral types for the companion star. 
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Figure 10. Equivalent width variations for H« and H/3, during 
Epoch B (JulOO minioutburst). 



speaking, this is a lower limit to the distance as we are ne- 
glecting any contribution by the accretion disc to the qui- 
escent optical flux. Although the measured (V-R) color is 
consistent with the calculated density, a spectral type de- 
termination of the companion star is essential to reflne this 
distance estimate. 



6.2 Relative contribution of the disc 

Following the color evolution of the system, we can compute 
the relative contribution of the disc to the optical flux in the 
V and R bands as a function of the assumed spectral type 
of the companion star. Assuming that the colour of the sec- 
ondary star corresponds to a normal dwarf of a fixed spectral 
type and that the colour of the disc remains constant during 
the main outburst (as we can see in Figure then we can 
relate the fluxes of the disc and star in the R and V bands: 
fi'disc— hxffl^.^^, fvstar^kxffl^j^^, where hoclO^*^~^''^"= and 
koclO'*'"^-'^''"'. It is straightforward to calculate the relative 
contribution of the disc (fv^ji^^/fy and fR^i^^/fn.) through: 

fv^hxf„, + {k-h)x fa,,^,. 

where and {r are the measured fluxes in the V and 
R bands. The unreddened colours ( V-R) corresponding to 
spectral types G5V to MOV are taken from Schmidt-Kaler 
(1982). We also assume E(B-V)=0.58 from Hynes et al. 1999. 
The results are shown in Table ^. 

During the first phase of the outburst, the optical emis- 
sion is almost completely dominated by the accretion disc. 
The relative contribution of the disc reaches a minimum 
when the system drops to the faint level between the June 
00 and July 00 minioutbursts. 
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6.3 Properties of the disc 

Using the same technique, we have fitted the unreddened 
B, V and R band fluxes to estimate the intrinsic colours 
of the disc. As we might expect, the fluxes during outburst 
are strongly correlated. We obtain {B-V)diac ~-0.07 and 
{V-R) disc ~-0.10. These colours agree well with irradiated 
model predictions, where most of the reprocessed energy 
is radiated in the UV (see e.g. van Paradijs & McClintock 
1995). 



The ratio of outburst X-ray to optical luminosity = 
Bq + 2.5 log Fx(/iJy)] agrees with the observed distribu- 
tion for LMXBs. Taking B=15. 9 (Chaty et al. 20 00) and 
Fx (2 - 12fcey) ~ 250 mCrab ( [Wood et al. 199E|) at th e 
outburst peak and assuming j4b=2.39 (Hynes et al. 199!:) 



we obtain ^=20.9 ±0.4, whereas the distribution peak of 
LMXBs gives ^=21.8 ±1 (see van Paradijs & McClintock 
1995). The moderately low ^ probably indicates a high bi- 
nary inclination since the X-ray source in J 1859+226 could 
be partially hidden by the accretion disc. This result is con- 
sistent with the evidence for two X-ray dips during the sec- 
ond outburst (~July 8) reported by Tomsick et al. 2000, 
and with the large amplitude (0.4 mag) of the quiescent el- 
lipsoidal modulation. 



6.4 The outburst mechanism 

It is largely accepted that optical emission in SXTs outburst 
is due to irradiation of the outer disc by X-rays. Within this 
context, the optical emission must be correlated with the 
X-rays. However, we see features which have no correspon- 
dence between the X-rays and optical, such as a secondary 
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maximum ~12 days after the peak of the X-rays outburst 
which is not visible in the optical band, and minioutbursts 
which are absent in the X-ray lightcurve. Such differences 
can not be explained in the context of a simple irradiated 
disc and an additional mechanism is necessary to explain 
them. Systems with c lear discrepancie s between X-rays and 
optical are J1655-40 ( |Esin et al. 2000[ ) and J1550-564 ( |jain 
et al. 230C). Both sources exhibit an exponentially declin- 



ing optical light curve, whereas the X-ray remains constant 
or increases slightly. To explain the outburst light curves 
in these systems it has been suggested the classical dwarf- 
nova type instability followed by an episode of enhanced 
mass transfer from the secondary or up scatter of the op- 
tical flux into the X-ray by the corona. Our photometric 
colors indicate that the optical emission is dominated by X- 
ray reprocessing on the disc although in J1859-I-226 clearly 
interplay with viscous heating and geometrical effects (such 
as disc height variations and shadowing effects) are playing 
an important role in the differences between the X-ray and 
optical lightcurves. 

The only SXTs where minioutbursts have been detected 
are XTE J0422-I-32 and GRS 1009-45, ahhough this prop- 
erty does not seem to be a peculiarity of these systems, but a 
selection effect, due to the difflculty in obtaining continuous 
monitoring from outburst to quiescence. The most promis- 
ing mechanism to produce minioutbursts is the X-ray echo 
model (Augusteijn et al. 1993, Hameury et al. 2000) , where 
they are interpreted as due to enhanced mass flow from the 
companion star on the outer disc which is responding, es- 
sentially linearly, to heating by X-rays from the primary. 
The process goes on continuously, triggered by echoing the 
initial X-ray outburst. 

The spectra taken during the main outburst show an 
emission feature at AA4540-4550A which is the Bowen Blend, 
but it is not present in the Jul 00 minioutburst. This behav- 
ior has also been observed in J0422-I-32 (Casares et al. 1995, 
Callanan et al. 1995). The Bowen blend is a combination of 
high excitation lines (mainly Cm, On and Nlll at AA4634- 
4642), produced by the fluoresencence resonance mechanism 
which initially requires seed photons of Hell Lya at A303.78. 
Under the hypothesis that the composition of the accreting 
material has not changed substantially since the onset of 
outburst, the absence of Bowen emission in the miniout- 
burst spectra must be ultimately related to the weakness 
of the X-ray photoionizing continuum which originates the 
whole cascade process. We do not have enough information 
in our data to make detailed calculations of this process. 
However, it seems that the reprocessed X-rays make a small 
contribution to the optical flux during minioutburst: in an 
irradiated disc, the average ratio of X-ray to optical lumi- 
nosity yields Lx(2-ll keV)/Lopt (300-700 nm)~500. We can 
estimate this rate during the highest minioutbursts (Jul 00) . 
Tomsick et al. 2000 report a 3.11x 10""erg cm"^ s"^ flux 
near the maximum of the minioutburst in the range 2.5-20 
keV, which was well above the expected quiescent X-ray flux 
level. Assuming the energy spectrum is described by a power 
law with a photon index 2.03 as they reported, we can esti- 
mate the flux in the range 2-11 keV. Taking a B magnitude 
of _B'~19.2, the estimated X-ray to optical ratio is Lx(2- 
11 keV)/Lopt (300-700 nm)~4. This means that the X-ray 



spectrum during minioutburst must be much harder than 
during primary outburst or, alternatively, that the optical 
flux during minioutburst is instead dominated by the in- 
trinsic luminosity of the disc. During the primary outburst, 
the intensity and width of the Balmer lines is notoriously 
smaller than in the minioutburst spectra as expected be- 
cause the continuum is 2 mag brighter. The change in line 
size is probably reflecting a shrinking in disc size. Alterna- 
tively the emission line cores might be filled in with narrow 
emissions during the main outburst. 



Another possibility could be that minioutbursts are pro- 
duced by a temporal enhancement of vi scositv in the quies - 
cent disc just after the main outburst ( Osaki et al. 1997 ). 
However, there is no easy explanation as to why the viscosity 
can vary in such a way. 
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